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Abstract. We give a brief overview and progress report on our lattice QCD calculation of neutral B mixing hadronic matrix 
elements needed for Standard Model and Beyond the Standard Model physics. Reference [1] contains more details and results. 
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1. MOTIVATION AND CONNECTION TO EXPERIMENT 



Neutral-B mixing is suppressed in the Standard Model, making it a process sensitive to new physics. The effective 
qq ' hamiltonian that describes neutral-B mixing J^g- = £? =1 Ci&i contains five operators whose matrix elements elements 

\0 '. (B^Wii^B^) must be calculated with lattice QCD, 
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<N ; 1 = (^y^Lq 01 ) (bPy^Lql 3 ), G 4 = (b a Lq a ) (PRql 3 ), 

02 = {b a Lq a ){PLqP) 7 <ff 5 = (b a Lql } )(bPRq a ), (1) 

ffi = (b a LqP) {PLq a ), 



where a and /3 are color indices. The first three operators are needed to describe Standard Model processes. The 
remaining two, G 4 and appear in Beyond the Standard Model physics. 
In the Standard Model, the meson mass difference can be expressed as 
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AM, = Mm b V t ;\ 2 {Sy^)\B° q ). (2) 

Measurements of AM q have sub-percent errors [2, 3], so our ability to to constrain the CKM matrix contribution V t bV* q \ 
is limited by how precisely we know (B^\Gi(ii)\B q ). Also of interest is the SU(3)-breaking ratio t,, defined by 
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E, is useful since some (lattice QCD) errors cancel in the ratio of matrix elements. Current lattice QCD calculations 
have errors of just under 7% on the square root of (B^\Gi (jj,)\B^) and just under 3% on | [4, 5, 6]. 
Including Beyond the Standard Model effects results in a generalized version of Eq. (2), 



am 9 = £q(m)(^I^(m)|b?}- (4) 

1=1 



Combined with lattice QCD results for all five matrix elements, Eq. (4) can be used to check that a Beyond the Standard 
Model prediction is consistent with experiment. The full set of matrix elements was calculated by Ref. [7] in 2001 
using the quenched approximation. Advances since then allow for an improved calculation. 
The lifetime difference, in the Standard Model, can be written schematically as [8, 9] 

Ar„ = fl t [G x B Xxj + G 3 B Xq ] cos 0, + 0(l/m b , a s ), (5) 

where fi, are bag parameters defined via (Z?|j | \B®) °< fg B; (ju ), the G,- are calculated perturbatively, and ty q is the 
CP-violating phase. In Eq. (5), Ar is dominated by the & \ contribution, resulting in a reduction of hadronic uncertainty 
in the ratio AT /AM. Nevertheless, a calculation of the ratio (B® | ^ 3 (il) |5°)/(B° | &\ (ll) on the lattice can be useful. 
In addition, Eq. (5) constrains the expected behavior of Ar versus [8, 10], and a calculation of certain combinations 
of the Standard Model operators should reduce the error in that relationship. 



2. THE CALCULATION IN LATTICE QCD AND A DISCUSSION OF ERRORS 

In these proceedings, we provide an overview of our error analysis; a more detailed description of the calculation can 
be found in Ref. [1]. The physics of a quantum field theory can be obtained from correlation functions which, when 
expressed in path integral form, can be evaluated numerically on a discretized volume of space-time (the lattice). These 
correlation functions can also be written as functions of matrix elements. There are a number of sources of error in a 
lattice QCD calculation. We review these and highlight improvements we have made over our previous calculation [5]. 

A lattice QCD calculation obviously introduces an error due to discretization. Guided by theory, we can extrapolate 
from a finite lattice spacing to the continuum and also quantify the discretization errors 1 . Our current lattice QCD 
calculation of the B-mixing matrix elements includes two lattice spacings, 0. 12 and 0.09 fm, included in our previous 
calculation as well as a third, smaller spacing of 0.06 fm. In addition, we plan to anchor our continuum extrapolation 
with a fourth spacing of 0.04 fm. In this way, we can gain better control over discretization errors. 

Numerical integration involves ensemble averages over gauge (vacuum) configurations. The number of gauge 
configurations available for each ensemble has increased by a factor of 2 - 4 for most of our data allowing for improved 
statistical errors. 

It is computationally expensive to calculate with up and down quarks at their physical mass, so we calculate with 
light quarks that are "too heavy". We then use chiral perturbation theory to extrapolate to the physical quark masses. 
Because the light-mass pions of XPTare the component affected by the finite volume in the simulations, the finite- 
volume error is folded into the chiral extrapolation. The range of "light" quark masses we use spans from the strange- 
quark mass m s down to 0.1m s , with one ensemble at 0.05m s . Such a range allows for good control over the chiral 
extrapolation and the decrease in statistical errors mentioned above improves these extrapolations. 

Lattice calculations also have errors from inputs: the physical scale of the lattice spacing and the bare quark masses. 
Each of these is determined by a separate lattice QCD calculation. Once determined, these inputs can be used for a 
suite of calculations of which B-mixing matrix elements are just one part. Details on the determination of these inputs 
can be found in Refs. [11, 12]. Finally, operators must be matched to the continuum and renormalized. For this we use 
one-loop perturbation theory resulting in a perturbation-theory truncation error. 

In the calculations discussed here, we generate numerical data for two-point and three-point correlation functions 
and fit them simultaneously for each meson (combination of heavy and light quarks). Once we have renormalized the 
matrix elements, we use SU(3), partially-quenched, staggered XPTfoi the extrapolation of each matrix element and 
ratios of matrix elements. 

Table 1 compares the error from our previous analysis to the error we expect to have for the calculation described 
here, based on our preliminary analysis. 



3. CONCLUSION AND SUMMARY 

We have a good start on a large-data-set lattice QCD calculation of the matrix elements that describe neutral-B mixing. 
Our calculation will cover the operators needed for both Standard Model and Beyond the Standard Model physics. For 



Light-quark discretization errors are folded into the chiral extrapolation; heavy-quark discretization errors are often quantified separately. 



the Standard Model matrix elements, we expect to halve the error on current, published calculations. In the Beyond 
the Standard Model case, this will be the first full-QCD (unquenched) calculation and the first update in ten years. 
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TABLE 1. Error budget on p Bs = ^ S\M E B in percent. Column "2008- 

09" is based on our calculation in Ref. [5]. "Expected" is our projection 
based on the preliminary calculation presented here. 



[!t] source 


2008-09 


expected 


Inputs: 






scale ij\ ) 


3.0 


1.1 


sea, valence quark masses (mis)tuning 


0.3 


0.3 


b-quark mass (mis)tuning 


1.1 


~ 1 


statistical 


2.7 


~ 1 


heavy-quark discretization 


2.0 


~ 1.2 


^PT + light quark discretization + finite volume 


0.7 


<0.7 


matching, renormalization (1-loop PT) 


~ 4 


~2.5 


total 


6.2% 


~ 3.4% 



